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Abstract 

^ We extend Kirman's model by introducing variable event time scale. The proposed flexi- 
ble time scale is equivalent to the variable trading activity observed in financial markets. 
Stochastic version of the extended Kirman's agent based model is compared to the non-linear 
stochastic models of long-range memory in financial markets. Agent based model providing 
matching macroscopic description serves as a microscopic reasoning of the earlier proposed 
stochastic model exhibiting power law statistics. 

Keywords: microfoundations, agent based models, stochastic models, financial markets, 
C/3 long-range memory 

d 

: i-H 
I 

i ^ 1. Introduction 

CN Computational modeling of complex systems pQ has become a rapidly developing method 

^ of doing science [2] . This approach is indispensable when microscopic nature of interactions is 
unambiguously defined. Nevertheless in the prevailing majority of social systems microscopic 

\Q 
(N 



nature of interactions has to be deduced from the macroscopic behavior of the whole system. 

There are two major types of financial market models - microscopic, most usually agent 
based, and macroscopic, usually based on stochastic calculus as most common choice in 
phenomenological description of complex systems. Agent based models (ex. [SHU [5]) come 
from the basic ideas about behavior inside the system. Resulting models can be rather 
.£h simple [5] or not so [3]. No matter the complexity of description there are still no agent 
based models, which would be able to sufficiently mimic empirical data [3J. In contrast, 
stochastic models (ex. [6J [7]) are being derived directly from the empirical data, thus 
reproducing observed statistical features, so-called stylized facts [HI [9]. The drawback of 
stochastic models is the lack of direct insights into the microscopic nature of replicated 
dynamics. Bridging between these two very different approaches one could propose model 
successful in both - theory and practice. 
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Top-down approach, starting from stochastic and moving towards agent based models, 
seems to be a very formidable task, as macro-behavior of complex system can not be un- 
derstood as simple superposition of varying micro-behaviors. In case of sophisticated agent 
based models [3] bottom-up approach provides too many opportunities. But there is a wide 
selection of rather simple agent based models (ex. [5]), whose stochastic treatment can be 
directly obtained from the microscopic description [TU] . 

In this contribution we consider the opportunity to generalize Kirman's ant colony model 
[5] with the intention to modify microscopic approach to the financial market pU] reproduc- 
ing main stylized facts of this complex system. In Section [2] we will briefly introduce Kirman's 
ant colony model, propose its generalization and derive stochastic model for population dy- 
namics. Next, in Section [3j we will introduce return, modulating return and propose single 
stochastic model for the modulating return. In Section [4] we compare generalized Kirman's 
model with stochastic description of financial markets analyzed in [TJ [TT] . Conclusion and 
discussion of future prospects will be given in Section [5j 

2. Kirman's model and its stochastic treatment 

In his seminal paper, [5], Kirman noticed that entomologists and economists observe 
similar patterns in rather different systems. 

Kirman credits entomologists Deneubourg and Pastels to be the first ones to observe 
ant colony with two identical food sources available (for the original papers see references 
in [5], while relevant research is also available in more recent paper [12]). They observed 
that in such case majority of ants still tend to use only a single food source at any given 
time. Though the other food source is not completely neglected as switches to the previously 
overlooked food source were observed despite the fact that food sources remain identical. 

Interestingly enough human crowd behavior tends to be quite similar, at least in sta- 
tistical sense. There are observations that majority of people tend to choose more popular 
product, than less popular, despite both being of a similar quality. The article [S] also cites 
numerous works, which speculate that herding behavior might be related to the endogenous 
fluctuations of asset price. 

Taking discussed empirical observations into account Kirman proposed a Markovian 
chain with the following transition probabilities: 



where N is a total number of agents (i.e. ants, traders or etc.) in the system (i.e. colony, 
market or etc.), X is a number of agents choosing one option (i.e. food source, trading 
strategy), while N—X use another option, At is a very short time step, for which p(X — » X+ 
1) +p{X — y X — 1) < 1 holds. In the above transition probabilities h terms describe herding 
behavior, as intrinsic property of agents themselves. While Oi terms describe individual 
transitions, which are assumed to be motivated by the attractiveness of the other option. 

Kirman also gathered empirical evidence that the aforementioned behavioral patterns 
are observed in ant colonies with very different populations - from tens to millions of ants. 




(N -X)(a 1 + hX)At, 
X[a 2 + h(N-X)]At, 
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This is most probably true in socio-economic scenarios. Thus it would be convenient to 
obtain a model description in the continuous limit, i.e. unrelated to N. One can assume 
that N approaches infinity securing continuity of x = in such case transition rates are 
expressed as follows: 



7T+(X) 



7T (X) 
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Evidently the above are related to the original transition probabilities, Eqs. ([T]) and (J2|, as 
p(X->X±l) = N 2 n ± (x)At. 

Using one-step process formalism [13] one can compactly express Master equation via 
one-step operators, E (increment) and E _1 (decrement): 



d t u(x,t) = N 2 {(E - l)[n-(x)uj(x,t)} + (E^ 1 - 1)[tv + (x)u(x, t)]} 



(5) 



where co(x,t) is time-dependent probability density function. As one-step operators act on 
continuous functions one can expand them using Taylor series up to the second order terms. 
After doing so one would obtain Fokker-Plank equation: 



d t uj(x,t) 

where in the limit of large N, 



-d x [A(x)co(x,t)} + ~d 2 x [B(x)cu(x,t)}, 



A(x) = N(7i + (x) — ir (x)) = cr 1 (l — x) — a 2 x 
B(x) 



tt + (x) + 7T (x) = 2hx(l — x) + ^(1 — x) + ^x m 2hx(l — x) 



This Fokker-Planck equation corresponds to the stochastic differential equation 



dx = [si(l - x) - e 2 x] dt s + J2x(l - x)dW, 



(6) 



(7) 



(9) 



where we have introduced dimensionless time scale, t s = ht, and accordingly scaled the 
original model's parameters, £j = y. In the above W s is appropriately scaled Wiener 
process. Eq. ([9]), in non-dimesionless form, was originally derived in slightly different 
manner by Alfarano et al. in [TO] . 



3. Stochastic Kirman's model for return in the financial markets 

As Kirman's model describes two state dynamics one must define two types of agents 
acting inside the market in order to relate Kirman's model to financial markets. Currently, 
the most common choice is assuming that agents can be either fundamentalists or noise 
traders [3]. 
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Fundamentalists are assumed to be long term investors who have fundamental knowledge, 
which is quantified as fundamental price, Pf(t), of the traded stock. Thus their excess 
demand, Df(t), is shaped by their long term expectations [TO] : 



D f (t) = N f (t) In 



m 
p( t ) 



(10) 



where Nf(t) is a number of fundamentalists inside the market and P(t) is a current market 
price. Being long term investors fundamentalists assume that P(t) will converge towards 
Pf(t) at least in a long run. Therefore if Pf(t) > P{t), fundamentalists will expect that P(t) 
will grow in future and consequently they will buy the stock (Df(t) > 0). In the opposite 
case, Pf{t) < P(t), they will expect decrease of P(t) and for this reason they will sell the 
stock (D f (t) < 0). 

Noise traders on the other hand are short term investors who estimate future price based 
on its recent movements. As there is a wide selection of technical trading strategies, one can 
simply assume that average noise traders demand is based on the mood, £(£), [10]: 



D c {t) = -r N c (t)£(t), 



(11) 



where D c (t) is a total excess demand of noise trader group, N c (t) is a number of noise traders 
inside the market and ro can be seen as a relative noise trader impact factor. 

Price and, later, return can be introduced into the model by applying Walrassian scenario. 
One can assume that trading in the market occurs trough the market maker, who sets a 
fair price. As fair price should stabilize market, sum of all groups' excess demands must be 
equal to zero: 



D f (t)+D c (t) 



N f (t)\n^~rM)U1) 



P(t) 



0. 



P(t) = P f (t)exp 



(12) 
(13) 



where without loosing generality one can assume that fundamental price remains constant, 
Pf(t) = Pf. Therefore return in selected time window T is given by: 



r(t) = In 



m 

P(t - T) 



r 



x(t) 



x(t - T) 



(14) 



where we have set that 



N c (t) 



x and 



l-x(ty y/ l-x(t-T) 

Nf ^ — 1 — x. Alfarano et al. simplified the above 



N M " ^ N 

by assuming that x(t) is slower process than £(t) [10J, obtaining adiabatic approximation of 
return: 

where ((t) — £(t) — £(t — T). If ((t) is modeled using spin-noise model, as in [TU], then 
can be seen as an absolute return. 



l-x(t) 
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Besides directly linking Kirman dynamics to the financial market scenario, we see a 
possibility to extend Kirman's model by using financial markets as an inspiration. It is 
known that trading activity (i.e. event rate) in financial markets is non-constant, while the 
original Kirman's transition probabilities, Eqs. ([TJ and ([2]), assume that agents meet (i.e. 
events occur) at a constant rate. As derivation of stochastic differential equation does not 
depend on the explicit form of n ± (x) (see previous section), we can easily introduce our 
extension directly to the Eq. (|9]): 



dx 



Ei(l -x) - 



e 2 x 
t(x) 



dt s + 



\ 



2x(l - x) 



-dW x 



T(X) 



(16) 



where t(x) adjusts time scale, of microscopic events, according to the current macroscopic 
system state, x. Note that E\ is not divided by r(x) - it is a consequence of assumed 
rationality of fundamentalists, namely their individual behavior should not be influenced by 
the current trades. 

Using Ito formula for variable substitution, which can be found in |14j . one can obtain 
stochastic differential equation for absolute return, y — • 



l-x{t) • 



dy 



£2 



t(v) 



[1 + y)dt s + 



I2y_ 

r (y) 



y)dW s 



(17) 



Note that absolute return, y in the above, serves as a measure of volatility in the financial 
markets. It is known that volatility has long-range memory and correlates with trading 
activity and has probability density function with power law tail [8]. 



4. Comparison with selected stochastic models 

In this contribution we study the case of r{y) = y~ a . This selection might be backed by 
the fact that trading activity has positive correlation with volatility. In such case obtained 



stochastic differential equation, Eq. (17), in the limit of y 3> 1 is very similar to the 
stochastic models discussed in [7] [TT]. 

In series of papers, see [TT] for the most recent discussion and relevant past references, 
stochastic differential equation, 



dy 



A N 



V 



y 



2r, ~ l dt + y v dW, 



was introduced as a class of stochastic differential equations providing solutions with power 



law statistics. It was shown that spectral density of stochastic variable defined by Eq. (18) 
is a power law: 

W-js- " = 1 + (19) 

Stationary probability density function p(y) of the aforementioned variable y is a power law 
as well 

p(y) ~ y~ x - (20) 
5 




Figure 1: Numerically calculated PDF (a) and power spectral density (b) of y denned by Eq. (171 in three 
distinct cases, a — (red squares), a — 1 (blue circles) and a — 2 (magenta triangles). Other model 
parameters were set as follows: E\ 
modelic results: (a) A 
magenta triangles). 



= 0, £2 = 2 — a. Solid curves are analytical fits, Eqs. ([24]) and ( |25| , for 
3 (red squares, blue circles, magenta triangles), (b) /3 — 1 (red squares, blue circles, 



Note that one has to introduce diffusion restriction, at least from the side of small values, 



of y in Eq. (18) 



such case Eq. (17) becomes: 



In the limit of large y, y 1, we can consider only the highest powers in Eq. (17). In 

3: 

dy = (2- e 2 )y 2+a dt s + J2y*+*dW s . (21) 



Direct comparison of Eqs. (18) and (]2l|) gives: 



n = 3 +" 
'/ 2 ' 

A = e 2 + a + 1. 



(22) 
(23) 



Consequently we expect stochastic variable y defined by Eq. (17) to have power law sta- 
tionary probability density function, 



— £2— a— 1 



and power law spectral density, 

S(f) 



p(y) ~ y 



j p 1 + a 



(24) 



(25) 



By using explicit form of spectral density, Eq. (25), we have reproduced 1/ f noise in three 
distinct cases of Eq. (17) (see Fig. [TJ. 

If we linearize drift function of Eq. ( JT7[ ) with the respect to return (i.e. set €2 = 2), we 
obtain stochastic differential equation (once again in the limit y ^> 1), 



dy = e iy dt s + ^/2y 3 + a dW s . 
6 



(26) 
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Figure 2: Numerically calculated PDF (a) and power spectral density (b) of y defined by Eq. (17) with 
linear drift function in three distinct cases, a = (red squares), a — 1 (blue circles) and a = 2 (magenta 
triangles). Other model parameters were set as follows: E\ — £2 = 2. Solid curves are analytical fits, Eqs. 
(28) and (29), for modelic results: (a) A = 3 (red squares), A = 4 (blue circles), A = 5 (magenta triangles), 
(b) j3 = 1 (red squares), j3 = 1.5 (blue circles), f3 = 1.66 (magenta triangles). 



which is similar to the generalized CEV process considered in [?]: 

dy = aydt + by v dW, 



(27) 



which is noted to be a special case of Eq. (18), when exponential restriction of diffusion is 
applied. Though comparison with this special case is important on its own as this equation 
generalizes some stochastic models used in risk management. Theoretical prediction of PDF 



and spectral density for y defined by Eq. (27), is given by [7]: 



S(f) 



p{y) 
1 

rs_/ — =■ 



y 

(3 



1 + 



a 
l+a> 



(28) 
(29) 



where we have used previously set relation between model parameters 77 and a, Eq. (22) 



As one can see in Fig. |2j these predictions are also correct for y defined by Eq. (17), in case 
of e 2 = 2. 

Financial market time series are also known to exhibit interesting scaling behavior, 
namely they are known to be multifractal [HI [15]. Thus a succesful model of financial 
market should generate multifractal time series. It is known that formaly constructed, via 
inverse Fourier transform, time series are not multifractal despite having proper spectral 
density, namely S(f) ~ iff 13 , while the point process based models were shown to exhibit 



both power law spectral density and multifractality [TB]. As Eq. (18) was primarily derived 
using the very same point process model we expect its and therefore Eq. (17) solutions also 
to exhibit multifractality. 

As one can see in Fig. [3] time series obtained by solving Eq. (17) are indeed multifrac- 
tal. Power law slopes, characterized by Hurst exponent, H(q), of generalized height-height 
correlation function, which is defined as 



F q (T) = (\v(t + T)-y(t)\*) 1/q , 
7 



(30) 
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Figure 3: Generalized height-height correlation functions, F q (r)^ versus lag, r, for the y time series obtained 
by solving Eq. (17l in three distinct cases, a — (a), a — 1 (b), a = 2 (c). Subfigure (d) is a plot of 
Hurst exponent, H(q), versus its order, g, - red curve corresponds to a = (or subfigure (a)), blue curve to 
a = 1 (or (b)) and magenta curve to a — 2 (or (c)). Other model parameters were set as follows: E\ = 0, 
£2 = 2 — a. In all three cases we have used single realization of million points to obtain F q (r). Nine orders 
of correlation, q = 1, 2, . . . , 8, 16, are ploted in each of subfigures. 



S 



where angular brackets, (...), denote time average, differ for different orders of correlation, 
q, in all three different cases of a. Note that for all orders of correlation Hurst exponent 
values are smaller than 0.5. This result is expected as mutlifractal time series exhibiting 
pink noise, S(f) ~l// /3 ,0</3<2, are known to have Hurst exponent smaller than 0.5 

usi- 



5. Conclusions and future work 

In this contribution we have started from a simple agent based model and obtained 



stochastic model for absolute return, Eq. (17). We found that obtained macroscopic model 



of absolute return is similar, in the limit of large values, y 1, to the empirically de- 
rived stochastic models considered in jTJ [11] . We have backed discussed analytical insights 
numerically - analytical predictions obtained in [7J [TTJ fit modelic results very well. 



As considered stochastic model, Eq. (17), was obtained from Kirman's agent based 



model, Kirman dynamics can be seen as a microscopic explanation of some very general 
non-linear stochastic models and, thus, statistical properties, namely long-range memory 
and fat tails, observed in actual financial markets. The exponents of power law statistics, A 
and (3, can be adjusted by introducing feedback between the macroscopic system state, y, 
and the rate of microscopic events l/r(y), where r{y) = y~ a . 

Basing ourselves on this simple model we can further work on replication of more so- 
phisticated statistical properties of return, namely fractured spectral density, and modeling 
of trading activity. We also think that it is possible to apply Kirman model more broadly 
- ideas underlying such birth-death processes are very well spread over different fields (ex. 
marketing [17] . social dynamic modeling [T5]). 
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